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1 Introduction 



Whilst considerable progress on multi-loop diagrams was made between the AI-HENP 
92 |T| and AI-HENP 93 |0] workshops, there is still a pressing need for methods that 
exploit hard-won analytical results by efficient techniques of numerical approximation. 
Progress in this direction was made, independently, in and [^], whose authors combine 
in this present work, which derives new analytical two-loop results for heavy-quark cur- 
rent correlators and uses them to refine the phenomenological extraction of the gluon 
condensate and to assess the numerical methods of 0, 

In Section 2, we obtain analytically, to two loops, the coefficient functions of the gluon 
condensate, (G^), in the correlators of heavy-quark vector, axial, scalar and pseudoscalar 
currents, for all values of z = q'^/Am^. This work is therefore an extension of that in ^, 
where the ffist 7 moments of these coefficient functions were analytically computed, by 
reducing the problem to the evaluation of vacuum scalar two-loop diagrams with one 
massless and two massive lines and arbitrary integer indices. Here we apply the more 
general strategy of 0, |^, reducing the problem to evaluation of two-point diagrams, of 
the type needed for the two- loop photon propagator [§] , and using the programs developed 
for to obtain explicit analytical expressions for such diagrams, in d = 4 — 2e spacetime 
dimensions. As input to these programs, we use the intermediate results of [|, H, obtained 
after the calculation of traces and the action of the projector for the operator G^. 

Renormalization of the resulting bare two-loop terms amounts to no more than mass 
renormalization of the d-dimensional one-loop contributions, which we perform in the on- 
shell scheme, expressing our results in terms of the pole-mass, m. Doing this, we discover 
that a term was missed in [§, when performing the limit e — in the MS scheme, 
with the effect that the results given in |^, ^ do not hold in the MS scheme (nor in 
any scheme that respects the Ward identity relating axial and pseudoscalar correlators). 
After correcting this inconsistency, we find agreement between our new results and the 
partial results of p, ^ for the ffist 7 moments. In addition to studying this z —>■ limit 
in Section 3, we also compare our results with those obtained in 0, as 2r — oo, and 
in 



rO|, as 2; — > 1, finding agreement with the leading term, in each case. 



In Section 4, we assess numerical approximations, developed in 0, |^, ^, [TT], |T2|, by 



comparing predictions, made on the basis of previously available input |^, |T^, with 
exact new two-loop results. The outcome is most satisfactory: with 7 moments as in- 
put, the methods achieve very high accuracy; with 4 moments, they fare almost as well; 
most remarkably, we shall show that a single moment, combined with only the leading 
behaviours as 2; — >■ 1 and z — >■ —00, gives the two-loop terms in the next 10 moments to 
better than 1% accuracy. 

Section 5 gives our conclusions, concerning existing phenomenology and future calcu- 
lations, both analytical and numerical. One significant motivation for studying numerical 
approaches with limited analytical input concerns the feasibility of approximating the 
three-loop photon propagator from a few of its moments, obtainable by the methods 
of IjT^, and then using it for applications such as the muon anomaly []T1|. We suggest that 



such an approach to ambitious calculations in both QED and QCD is indeed feasible. 
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2 Exact results 



We evaluate contributions to the correlators, i / e"'^'(T( J(x) J(0)))dx, of the vector (V) 
current = ip'j^ip, the axial (A) current = V'TajTsV'! the scalar (S) current = 
2'mip'il', and the pseudoscalar (P) current = 2mipij5ip = d'^J^, where ip is a. heavy- 
quark field of mass m. We choose to work in the on-shell (OS) scheme, where the results 
are simplest. Since the currents have no anomalous dimensions, one may translate the 
results to any other scheme (e.g. the MS scheme) by making a one-loop transformation 
from the pole mass, m, to the renormalized mass of that scheme (e.g. by using m = 
m(/i)(l + (ln(yU^/m^) + |)as/7r + O^a"^)), where m(/i) is the MS mass, at scale fi). 

We denote the correlators of J^'^ by 11^'^ (g^) and decompose the tensor structure of 
the vector and axial correlators as follows: 

i/e^'^^(T(jJ(a:)jJ(0)))dx = n^(g^)(g,g. - g V) , (1) 



where 11^(0) = —A{m%l)%l)) enters (Q) as an equal-time commutator |15, 16|. Finally, we 



define dimensionless coefficients of the non-perturbative gluon condensate by writing the 
dimension-4 contributions to the correlators as 

KM') = ^^"^(1?:,^°^^ (C'iz) + OiaD) ; C\z) = C({z) + ^Ci{z) , (3) 



^°P''' ' (2m)"./ V ^ ' ' ^ ' ^ ^ ' Ti 

in the channels J = V, A, S, P, with exponents rij = 4, 2, 0, 0, respectively, making C'^(z) 
a dimensionless function oi z = q^/Am?, with one- and two-loop contributions C(2{z)- 
Since the currents are not renormalized, and UsG"^ is not renormalized at one loop, the 
renormalization scale of as in the two-loop term is irrelevant here. It should, however, be 
taken as 0(m), to suppress large logarithms at three- loop order [Q. 

Using Reduce for the trace calculations of [^, Q and for the recursion of the 
resulting scalar integrals of Q], we obtained the following on-shell results in the 4 channels: 

A^z\l- zfC^{z) = -?> + Az- Az^ + ?>{!- 2z)G{z) + —P^{z) , (4) 

vr 

Sz{l-z)C^{z) = l-2z-G{z) + —P^{z), (5) 

TT 

lQz{l- z)C'^{z) = -l-2z + {l-Az)G{z) + —P^{z), (6) 

vr 

AS{l-zfC^{z) = 7-10z + 3(3-4z)G(^) + — PP(;2), (7) 

TT 

P\z) = P({z) + P:^{z)G{z) + Pi{z)G\z) + {l-z)Pi{z)H{z), (8) 
with polynomials, Pf{z), given in Table 1. As d —>■ 4, two basic integrals are encountered: 

2v a/1 - 1/z - 1 
G{z) = ^Iny; y = \ \ , (9) 

H{z) = -^(Li2(2/2)-Li2(2/) + (21n(l + y) + ln(l-y)-flny)lny), (10) 
with dilogarithms Li2(?/^) = Z]n>o l/^"/""-^! P = 1;2, giving Clausen's integral if z G [0, 1]. 
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3 Limiting cases 



For 2; — > 0, one may approximate the coefficients by truncation of their Taylor series 



C-\z) = j:(< + ^ci)z-; bi^ci/ai, (11) 

n=0 ^ 



with one-loop moments given by 

V _ 2n + 2 (4)„ ^__li3)n s _ 3n + 4 (2)^ p _ _?7^ (2)^ 

15 (!)„' 3(f).' 12 (§)„' ^" 12 (!)„' ^''^ 

where (a). = r(a + n)/r{a). The two-loop corrections, b:^, are easily obtained from 
Table 1, using the expansions |1T6[| 



00 oo/n-| ""^X 

G{z) = E , (1 - z)G\z) = E > ^(^) = E E ^ - E I ^n-" , 

n=0 n=0 "I" n=0 \r=0 ^' "r -L r=l ^' / 

(13) 

with (7„ = n!/(3/2)„,. Results for n < 10 are given in Table 2, which may be extended up 
to = 200 in a few minutes of CPUtime. Since = 0, 64 is undefined. We find that 
C4 = —3266/4725. Two features of Table 2 are notable: the relation b^ = 6^ ensures the 
absence of a singularity in (^ as — > 0; the scheme- independent value feg = 6q = 11/4 
gives, via a Ward identity |jT5|, |16| , the new two- loop term in the heavy-quark expansion 

W^) = = fl + + Oial)) + 0(l/m^) , (14) 



127r V 4 vr 

whose one-loop term was used in |TP|. In the MS scheme, at fx = m, one obtains two- loop 

corrections Tk^ = bf^ — |(2?2 + rij). Comparing Table 2 with the corresponding MS tables 

in 0, we find that the latter give (6^ — 2), instead of 6^, for the reasons given in the 
introduction. 

Next we consider the limit z —00, which yields the asymptotic corrections 

_^^4 + 21„^-,j. (15) 




To relate these results to previous work, we use the methods of fl^ . The vector and axial 
light-quark results of 0, combined with our heavy-quark result (|14D, give 

IZTT \ b H J \ 6 H J IZTT V TT 



A , 9, {oisG"^) ( 7as\ , ^ ,^ ( 7a;s\ (asC^) / 34as 



The (pseudo) scalar results of |^ similarly agree with (|T^JT^). Note, in (|TB|), the change in 
sign between one-loop light- and heavy-quark terms. (The authors of H regret that a 



one-loop sign error in their previous work obscured the need to use [19 



Finally, for the threshold behaviour in the vector channel, as z — > 1, we obtain 



r^(^\ — 2304^' I 768 ^' 6912^' , 72^' ^^^-^ ' \ Cl \ I 

^2l^J- (1 _ + (1 _ ^)5/2 (1-^)2+ (1-^)3/2 +^1^(1-^)3/2^' 

whose first term agrees with |T0[. (An input used in is ruled out by the second term.) 
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4 Approximations 



We now assess previous numerical methods ||, H, 0] , by testing the accuracy to which 
they predict new features of the exact result for Cj(z). As input for the methods we take 
the following 'old' data: the leading asymptotic behaviour, from f^; the leading term in 
the threshold expansion, from the radiative corrections {b^\n < 7} to the first 7 
moments. (This is the input used in |^, corrected for errors in @, |^, |^.) 

First we consider the approximate spectral ansatz of [§, which is identical to using 
the following (exact) one-loop and (approximate) two-loop threshold expansions: 

N+3 

C^{z) = -lF2{z) + l^F,{z), Cj(z)^ ^/,F,(;.); Ffc(^) = 2^1(8, 1; ^; z) , (19) 

k=l 

with hypergeometric basis functions, Fk{z), giving polynomials in 1/(1 — z), for k = 1,3, 5, 
and terms involving G{z) or ln(l — z), otherwise. The {N + 3) two-loop coefficients, fk, 
are determined by the first N moments and by 3 further constraints [Q]: 

1977r2 ^+3 10 

^^ = -^' T.ik-m = 0, Y.k{k-l)U = -, (20) 

^-^U^ k=2 



giving the leading threshold singularity, derived from |T^, and z^C^(2;) — > —5/12, derived 



from IP], as z — >■ — 00. We remark that ansatz (|T9|) of P| cannot reproduce the form of 
the 4th term in the true threshold expansion (p!8|). We assess it by comparing the values 
of /2,3, required by the input, with those required by the second and third terms in 



namely /2 = 65/144 and /s = — 4137r^/3456, and also by comparing the output for 
{6^1 n = N . . . 10} and B{z) = C^(2;)/C^(z) with exact results. Table 3 shows the very 
high accuracy achieved with 7 moments, in all but the /2,3 tests. The results using only 
4 moments are almost as good, except for 2; ~ 3. Remarkably, the = 1 column shows 
that just 3 input numbers, namely , fi and B{—oo), give 10 additional moments to 
better than 1% accuracy, which is most encouraging for three-loop applications. 

For the last 3 columns of Table 3, we mapped [|12| and Fade- approximated p|, H, |12| 



z{l-zrcnz) + --^ = ^; ^=^^, (21) 



with the z-plane, cut along z G [1, 00], mapped to the unit disk, |ci;| < 1. By construction, 
D{uj) is finite at = 1 (i.e. z = 1) and diverges only logarithmically as — » — 1 (i.e. z 
—00). Its value at = (i.e. z = 0) is —fi- With = 1,4, 7 moments as input, we 
computed the [0/1], [2/2], [3/4] Fade approximants to D{uj), which were expanded to give 
the further moments of Table 3. From the approximations to -D(l) and -D'(l), we obtained 
the second and third terms in the Fade-approximated threshold expansion and compared 
them with the exact ones, in (0). Finally, we computed Fade-approximated values of 
B{z) = C2{z)/CY{z), using [|l2] u = exp (^2i arccos 2""'^/^) in the cases with z G [1, 00]. 

Table 3 clearly demonstrates that both (|19D and ( ^1]) are highly effective predictors of 
additional moments, even with very limited input. More input is needed to achieve high 
accuracy on the cut, since neither method performs well in the /2,3 tests. With A^ = 10 
moments, the Fade method gives the modulus of Cj(2;) to 0.05% and the phase to 0.04°, 
on the entire cut, whilst the figures for ansatz (|T9|) are 0.5% and 0.2°, respectively. 
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5 Summary and conclusions 



Our main analjd;ical results are in (|4HlO|) and Table 1, which give the two-loop coefficient 
functions of the gluon condensate in the correlators of heavy-quark currents. As suggested 
in H , computer algebra in d dimensions proves to be a most efficient way of obtaining 
new 4-dimensional results, whose analytical complexity enters only at the final stage, in 
this case via the dilogarithms of (p!0|). 

By combining and refining the REDUCE programs developed for 0, ^, we obtained 
a complete analytical result in the vector channel in about 2 hours of CPUtime on a 
486DX2-50 PC, whose timing in the standard Reduce test is 3 s. In comparison, it took 
20 hours to obtain merely the ffist 7 moments in 0, using a machine with a benchmark 
of 11 s. Our improvement in generality and speed entailed considerable programming 
effort to produce efficient code, of the type used in |Q, in order to process the many 
scalar integrals in the 0.7 Mbyte of output from the programs used in [^]. This involved 
systematic implementation of recurrence relations for 5F4 hypergeometric functions, which 
was exhaustively checked by independent Form code, written for 0]. 

We have found full agreement with previous results in three limiting cases: z —>■ 0, 
z ^ —00, and z — > 1, after correcting errors in [§, ^, P]. Using a Ward identity, we have 
obtained the new two-loop term in the heavy-quark expansion ([T^). 

Combining our new results with the procedure used in 0, to estimate the gluon 
condensate from empirical data ||2^, we arrive at {{ois/T^)G^^j^G^^) ~ 0.021 GeV^, as 



compared with 0.025 GeV^ in [Q. Our two-loop extraction thus gives a value that is still 
about twice as large as the effective one-loop value 

We have critically assessed two methods of numerical estimation: the approximate 
spectral ansatz ([T9|); and Fade approximation 0, |12[ after the mapping (^). Table 3 



shows that each performs rather well, even with limited input. Despite its analytical 
simplicity (rational approximation to D{uj) in (0)) the Fade method is as good as ( ]T9| ) 
for moment prediction and better for achieving high accuracy on the cut. We recommend 
both methods for future applications, with agreement between them providing a useful 
cross check. Note that the Fade method is more general, as it does not rely on the 
specific form of a lower- loop result. Most importantly, one should use information from 
all 3 limiting cases, in each method. We find that 4-figure accuracy can be obtained on the 



cut using only 10 moments in the Fade method of (^), whilst 17 moments are needed to 
achieve comparable accuracy using only the methods of [Q, without the limiting values 
as z ^ 1 and z —00. 

In conclusion, we believe that our analytical progress and accuracy of numerical ap- 



proximation offer real prospects of extending the two-loop approaches of 0, [TT|, |T2[ to three 
loops. Future work can be based on calculation of a few three- loop terms as z ^ —00 
and z ^ ||1^, |1^, for both QED and QCD two-point functions. In both theories, 
the leading three-loop term at threshold is already known |jTO|. We anticipate being 



able to obtain an accurate three-loop approximation for the photon propagator, for use 
in multi-loop QED [Q, and accurate estimates of three-loop corrections to QCD sum 



rules. Finally we note that mapped Fade methods are well suited to a very wide variety 
of standard-model applications, as indicated by their successful application to two-loop 
three-point functions P, [12|. 
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Table 1: The polynomial coefficients of ®, for use in 

pV(^) ^ -^(613 + 88282 - 22092^2 ^ 13248^3) 

p^(z) = y^(659 - 3774z + 30184^2 _ 55344^3 ^ 26496^^) 

P^(z) = -^(23 + 1456z - 41064^2 + 195024^^ - 350688^^ + 275328;^^ - 79488^^ 
p^[z) = 1^(473 - 3614^ + 6504^2 _ 33i2;z3) 

P^(z) = 216(1241 - 8390^ + 6624^2) 

P^(z) = -2Yg(1339 + 4128;z - 17656^^ + 13248^3) 

P^(z) = j^{A9 - 24042 + 32328^2 - 98112^3 ^ io7616z^ - 39744^^) 

P^(z) = ^z(467 - 2OOO2 + 1656^2) 

Pf(z) = -2i(77 + 2402z - 2208^2) 

P|(z) = ^(187 - 1636z + 5416^2 _ 4416^3) 

pS(^z) = -Y^(55 + 82^ - 8880^2 + 30624^3 _ 35040^^ + 13248;25) 

P|(z) = |z(121 - 6322 + 552^2) 

pP(z) = -^(3213 - 95422 + 6624^2) 

pP(z) = ^(2739 + 94282 - 261042^ + 1324823) 

pP(z) = -3^(99 - 114702 + 73008^2 - 1517762=^ + 1300802^ - 397442^) 
pP(z) = 1(54 - 12992 + 293622 - 16562^) 

Table 2: Values of = ci/a^ for J = V, A, S, P and n < 10 



n 

1 
2 
3 
4 
5 
6 
7 
8 
9 

10 



V 


A 


S 


P 


1469 


131 


11 


11 


162 


24 


4 


4 


135779 


3193 


2447 


131 


12960 


486 


504 


24 


1969 


72631 


2813 


307 


168 


9720 


450 


45 


546421 


124847 


160591 


97073 


42525 


15120 


21840 


15120 


661687433 


36591619 


400895 


undefined 


47628000 


4082400 


48384 


1800293669 


604286927 


1457 


28678327 


121080960 


62868960 


160 


1663200 


180657583657 


5566166237 


16727257 


1854828721 


11442150720 


544864320 


1698840 


113513400 


72346752463 


26439650563 


398248999 


72267337 


4341887550 


2451889440 


37837800 


4324320 


377046256819 


9437661186719 


321477183403 


5495760533 


21549939840 


833642409600 


28817268480 


316673280 


22322605695461 


98175636210793 


30518020922371 


159068122366537 


1220106888000 


8296726838400 


2597964969600 


8799558768000 


162570424650982 


38340293613343 


2191394981453 


28163524523 


8527479553593 


3111272564400 


178086308400 


1496523600 
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Table 3: Assessment of ansatz (|T^ and mapping (^Tj), using moments 



test 


exact 


ansatz 


ansatz 


ansatz 


mapped 


mapped 


mapped 


feature 


result 


AT ^ 


AT A 

iV = 4 


AT 1 

iV = 1 


AT ^ 


AT A 

iV = 4 


AT 1 

iV = 1 




9.06790 


input 


input 


input 


input 


input 


input 


bl 


10.4768 


input 


input 


10.4484 


input 


input 


10.4778 


b^ 


11.7202 


input 


input 


11.6675 


input 


input 


11.7281 


bs 


iZ.o4y4 


input 


input 


iZ. M Oi 


input 


input 


iz.oooo 


iN 
O4 


13.o928 


input 


13.8928 


13.8019 


input 


13.8927 


13.9203 


/iV 

»5 


i4.o0o0 


mput 


i4.o0o4 


14. r OzO 


mput 


1 /I fiRSO 
i4.o0oZ 


1 /I onfii^ 
i4.yuoo 


bj 


15 7888 


input 


15 7885 


15 6695 


input 


15.7881 


15 8369 




1 6 6625 


1 6 6625 


1 6 6620 


1 6 531 5 


16 6625 


16 6614 

XVJ . VJVJ Xt: 


1 6 7204 




1 7 4qfi4 


1 7 4qfi4 


1 7 4955 


1 7 3549 


1 7 4qfi4 


1 7 4q4fi 


1 7 5636 


hH 
"9 


1 8 2956 


18 2956 


18.2944 


18.1446 


1 8 2956 


18 2931 


18 3715 

X . 1 X(J 


"10 


1 9 0643 


1 9 0643 


19 0626 


18 9046 


1 9 0643 


19 0609 


1 9 1 485 


-Jl 


0.84388 


input 


input 


input 


input 


input 


input 


J2 


0.45139 


0.45868 


n rn /I on 

0.50439 


n r ^ 1 1 

0.57116 


0.47972 


n r 1 /I 

0.53134 


n /I n'7n cr 

0.40795 


P 

-13 


1.17944 


1.49727 


2.09688 


00 '70 

2.23783 


1.86522 


onon/^ 

2.39296 


1.62672 


5(-CX)) 


5.00000 


input 


input 


input 


input 


input 


input 


Bi-9) 


5.58641 


5.58639 


5.58495 


5.62888 


5.58733 


5.58208 


5.61983 


B{-3) 


6.28843 


6.28843 


6.28800 


6.33614 


6.28847 


6.28730 


6.31060 


B{-0.9) 


7.41565 


7.41565 


7.41561 


7.44826 


7.41565 


7.41557 


7.42198 


B{-0.3) 


8.29735 


8.29735 


8.29735 


8.31290 


8.29735 


8.29734 


8.29824 


5(0.3) 


10.3704 


10.3704 


10.3704 


10.3449 


10.3704 


10.3704 


10.3745 


5(0.9) 


25.6939 


25.6939 


25.6835 


25.4851 


25.6928 


25.6745 


25.8290 


|5(3)| 


1.81071 


1.80980 


1.91045 


2.47868 


1.80718 


2.00465 


2.21432 


|i?(9)| 


4.47078 


4.47366 


4.44784 


4.46713 


4.46202 


4.43571 


4.69021 
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